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ABSTRACT 



We present total and final-state resolved charge transfer (CT) rate coefficients for low-charge Ge, Se, Br, Kr, Rb, and Xe ions reacting 
with neutral hydrogen over the temperature range 10^-10* K. Each of these elements has been detected in ionized astrophysical 
nebulae, particularly planetary nebulae. CT rate coefficients are a key ingredient for the ionization equilibrium solutions needed to 
determine total elemental abundances from those of the observed ions. A multi-channel Landau Zener approach was used to compute 
rate coefficients for projectile ions with charges q = 2-5, and for singly-charged ions the Demkov approximation was utilized. Our 
results for five-times ionized species are lower limits, due to the incompleteness of level energies in the NIST database. In addition, 
we computed rate coefficients for charge transfer ionization reactions between the neutral species of the above six elements and 
ionized hydrogen. The resulting total and state-resolved CT rate coefficients are tabulated and available at the CDS. In tandem with 
our concurrent investigations of other important atomic processes in photoionized nebulae, this work will enable robust investigations 
of neutron-capture element abundances and nucleosynthesis via nebular spectroscopy. 

Key words, atomic data - atomic processes - Hll regions - nucleosynthesis, abundances - planetary nebulae: general - stars: evolution 
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The detection of n-capture elements (atomic number Z > 30) 
in ionized astrophysical nebulae has underscored the need for 
atomic data required to accurately solve for the ionization equi- 
librium of these species. Trans-iron element emission lines 



fe.e.. Zhane et al.l2005: Sharoee et al. 2007 


; Sterling et al.l2007t 


[Sterline & DinersteinI 20081 ISterlina et al. 


2009; Otsuka et all 



|2010, 2011) to H II regions (e.g., Baldwin et al 2000^ 
[Blum & Mc GregoJ l2008l and references therein), and th e 
interstellar medium of other galaxies (IVanzi et al.l l2008h . 
Determinations of their abundances have important implica- 
tions for the production and chemical ev olution of «-capture 
nuclides in th e Unive rse (e.g ., Busso et a l. 1999; Sneden et al. 



e.g.. muss 

l2008t ISterUng & DinersteinI |2008|) . Se, Kr, and Xe are the 
mos t wid ely detected n-capture element s in ionized nebulae 
re.g. | pine rstein 2001; Sharpee et al. 2007; Sterling & DinersteinI 
l2008t Isterling et al.i 12009) . but other trans-iron elements 
have also be en de t ected in nebular spectra, including Ge 
(Sterling et a T]l2002t ISt erling & Dinerstein 2003), Br, and Rb 



ionization corrections, provided that accurate data for atomic 
processes affecting the ionization equilibrium of the modeled 
elements are available. Since «-capture element emission lines 
have been detected primarily in photoionized nebulae, the most 
important atomic data required for ionization balance determina- 
tions are photoionization cross sections and rate coefficients for 
radiative recombination, dielectronic recombination, and charge 
transfer (CT). Unfortunately, these data are poorly if at all known 
for ions of trans-iron elements. 

This paper is the third in a series presenting theoretical 
atomic data determinations relevant to nebular ionization bal- 
ance solutions for low-charge n-capture element ions. In the first 
two (Sterling & Witthoeft 2011; Sterling 2011), we furnished 
multi-configuration Breit-Pauli distorted-wave photoionization 
cross sections and radiative and dielectronic recombination rate 
coefficients for Se and Kr ions. 

To address the need for CT data for low-charge ions of n- 
capture elements, we present calculations based on the Landau- 
Zener (LZ) and Demkov approximations. We consider reactions 
of the form 



(Pequignot & Bal uteaul ll994i; ,Zhang et al., ,2005i: iSharpee et al.1 
[2007; Steriing et al.l l2009l) . 

Due to the limited number of detectable n-capture element 
ions in individual nebulae, it is necessary to convert abundances 
of the observed ions to elemental abundances by correcting for 
unobserved ionization stages. Numerical simulations of ionized 
nebulae present a robust and accurate method of deriving such 



X«+ H- H ^ X*''""^ H- H"- -H AE, 



(1) 
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with q = 1-5 and X = Ge, Se, Br, Kr, Rb, and Xe. The forward 
reaction, CT recombination (in terms of the projectile X''^), is 
exothermic for these elements when q > 2, and the reverse re- 
action for multiply-charged systems is generally unimportant at 
the temperatures of photoionized nebulae (near 10'* K). On the 
other hand, of these elements only Kj has a higher ionization 
potential than H, and hence CT recombination is endoergic for 
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singly-ionized Ge, Se, Br, Rb, and Xe. We therefore provide rate 
coefficients for both CT recombination of singly-charged ions 
and CT ionization of neutral species. 

Our LZ calculations are multi-channel, in the sense that 
several exit channels for the X**"''^ system are simultane- 
ously considered. This avoids the problem of overestimating 
cross sections when adding two-channel cross sections of the 
same molecular symmetry (Butler & Dalgarno 1980b; Kingdo^ 
Il995h . For singly-ionized atoms, CT usually cannot be treated 
in the LZ approximation, and hence we utilize the Demkov ap- 
proximation to determine rate coefficients for those systems. 

LZ calculations are generally accurate to within a factor of 
three for systems (and final channels) with large rate coefficients 
(> IQ -^ cnP s-' HButler & DalgarnolfT980bl: iKingdon & Ferlandl 
[l996h . but can be less accurate when the rate coefficient is small 
(though in that case, CT often is much less important than ra- 
diative and dielectronic recombination). The quantum mechan- 
ical molecular-orbital close-coupling (QMOCC) approach (e.g. , 
iKimura & LanellT989t IZvgelman et al ] ll992l:IWangetal.ll2004 
is the most accurate theoretical method for determining CT rate 
coefficients at the temperatures of interest. However, given the 
time- and computationally-intensive nature of such calculations, 
we have chosen to first render LZ calculations for these systems. 
We will test the sensitivity of abundance determinations to un- 
certainties in the CT rate coefficients (assumed to be roughly a 
factor of three) through numerical simulations of ionized nebu- 
lae. Those systems which are most critical for ionization equilib- 
rium solutions are targets for further investigation with QMOCC 
methods. 

The structure of this paper is as follows. In Sect. |2] we dis- 
cuss the details and methodology of our calculations. We present 
the resulting total and final-state resolved CT rate coefficients in 
Sect. [3] and summarize our investigation in Sect. |4] 



2. Calculations and methodology 

2.1. Multi-channel Landau Zener calculations 

For projectiles with charge q - 2-5, we compute CT rate co- 
efficients in the multi-channel Lan dau-Zener (M CLZ) approxi- 
mation, based on the formalism of iButler & Da lgarno (1980b) 
and Jane V et al. (1983) with the low-energy modifications of 
Iciiang & Pritchai-d (.19791) . 

WTiile CT is a quasi-molecular problem, the basic details of 
the interaction can be understood from consideration of the ini- 
tial channel (in which both the projectile X*^ and target H are 
assumed to be in their ground states), the final channels (which 
can include ground and excited states of X**"''^), and the inter- 
action between them. As the atoms approach each other, if the 
ionization potential of X**"'*^ exceeds that of H (as is the case 
for q - 2-5 for all elements we consider), CT is exoergic and the 
diabatic potential curves can cross at an internuclear distance R^. 
The CT reaction is usually driven by a radial coupling that con- 
nects initial and final channels corresponding to the same molec- 
ular symmetry of XH''^. Adiabatic potential curves of molecular 
states with the same symmetry cannot cross, and R^ corresponds 
to the radius at which the adiabatic potentials most closely ap- 
proach (the "avoided crossing" distance). 

In the LZ approximation, charge exchange is assumed to be 
localized at the internuclear distance R^, with the probability of 
exchange given by 2p(l - p), where p = e"''' and 



n^[AU(Rdf 



hv[i(Hu-H22)l 



(2) 



In this equation, v is the relative radial velocity of the projectile 
and target, H[ i and H22 are the incoming and outgoing diabatic 
potentials, and AU (R^) is the separation of the adiabatic poten- 
tials at R^. The cross section at energy E is 



o-iE) = AnRipoil +A) \ (1 - e-'")e-'"x-^dx, 



(3) 



where A = [i/ii(oo) - Hu(R^)]/Ei ("Butler & Dalgarnolll980bh . 
The approach probability po is determined from the orbital and 
spin angular momenta of the initial and final states, which can 
form gif = (2Li + 1){2S i + 1)(2L/ + 1)(25/ + 1) molecular 
states, and the spin S and orbital angular mo mentum projec- 
tion quantum number A' of the molecular state (iHerzberd 19521: 
iDalgarno et aL.1990 ): 



po = (25 + 1)(2 - 6o,A')/8if. 



(4) 



The above describes a two-channel LZ interaction. In MCLZ 
calculations, several capture channels are considered simulta- 
neously, each with a unique value of R^, as the projectile ap- 
proaches the target. In this case, for fin al channels the pro ba- 
bility for capture into a state n is given by (I Jane v et al.lll983h 



Pn = piP2---P„(l - Pn)['^ +(Pn+\Pn+2---PN) + 
(Pn+\Pn+2 ■ ..pN-lf(l - Pn)^ + ■■■ 
+/'«+! (1 - Pn+lf + (1 - Pn+\fl 



(5) 



LZ calculations require estimates of the incoming and outgo- 
ing potentials, the avoided crossing distance R,^, and the energy 
separation of the adiabatic p otentials at R^. We use appr oxima- 
tions for these quantities from lButler & Dalgarnd(ll980bl) as fol- 
lows. In atomic unitsQ the incoming potential is 



Hn{R) = 



2R* 



+ Aqe 



-(0.&+()R 



(6) 



where is the polarizability of H, q the charge of the projec- 
tile, R the internuclear distance, and ^ is the exponent of a single 
orbital wavefunction (-1.0 for H targets). Typically set to 25, 
A was reduced in some cases to yield physically reasonable be- 
havior in the potential, which occasionally was spurious in the 
repulsive core (this generally only occurs when 7?x ^ 4.5-5 ao, 
in which case the rate coefficient is usually small at photoionized 
plasma temperatures). The outgoing potential is given by 



H22iR) = ^^-^ - AE, 



(7) 



where AE = IP(X'?+) - IP(H) - Eexc is the energy defect, IP is 
the ionization potential of the indicated system, and fiexc is the 
excitation energy of the final X'^"'^^ state. 

The avoided crossing distances R^ are determined by equat- 
ing Equations|6]and|7] In Fig.[Tl we illustrate the diabatic poten- 
tials and R^ values for the final channels of the case X^^ - Kr^^. 
Note that we ignore fine structure in our calculations. 

The final piece of necessary information, AU{R^), is perhaps 
the great est uncertainty in LZ cal culations, as the formula pro- 
vided bv Butler & Dalgarno! (1 1980 b) is an empirical fit to a small 
number of quantal calculations. The value of AU (R^) depends on 
whether the electron capture process is Type I (the H 1 5 electron 
is transferred without affecting the core electrons of the projec- 
tile) or Type II (electron capture is accompanied by the excitation 



' All quantities in this paper are in atomic units unless otherwise 
noted. 
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4 6 8 10 12 14 16 If 
Internuclear Distance (a.u.) 

Fig. 1. Incoming diabatic potential Hn for Kr-^^ + H (thick 
solid line), and diabatic potentials H22 for six exit channels of 
^ ^. corresponding to and ^S" molecular sym- 
metries (all listed states have an [Ar] 4s- 4p^ core). The intersec- 
tions of the incoming and outgoing potentials correspond to the 
avoided crossing distances for each exit channel. The poten- 
tial energies are shifted by taking Hii(oo) - 0. 



of a core electron). The formulae of iButler & Dalgamol (Il980bh 
for these two scenarios are: 

AUi = 27.2i;?^g-/3Rx (8) 
AUii = O.Se-^'^'^-AUi, (9) 

where /3 - 2IP(H) - 1.0 in atomic units. The potential for 
Type II reactions is reduced relative to their Type I counterparts 
due to the rearrangement of the core electrons. 

Once the quantities Hn, H22, Rx, and AU{Rx) are deter- 
mined, then the cross section for each final channel can be deter- 
mined from Equation [3] as illustrated for X''^ = Kr^^^ in Fig. |2l 
These are integrated over a Maxwellian velocity distribution to 
derive final-state resolved rate coefficients. 

2.2. Demkov calculations 

In the case of singly-charged ions (q = 1 in Equation [T]l, the 
LZ approximation cannot be applied since the final channels do 
not have Coulomb potentials (hence there are no avoided cross- 
ings). The reaction proceeds through radial coupling connecting 
molecular states with the same symmetry as in the LZ calcula- 
tions, but the coupling is not localized and can be represented 
with an exponential form. This type of interaction can be de- 
scribed by the Dem kov approximation dPemkovl 119641 : lOlsonI 
ll972l:ISwartzill994 . 

The reaction is assumed to occur at an interaction distance 
that we label R^ due to its analogous nature to the avoided cross- 
ing d istance in LZ calculations. /?x is given by equating AU{R) 
with (ISwartzlll994 



AH(R) = i/ii(X+ -H H) - H22(X + H+) - AE. 



(10) 



For the potential //ii(X -1- H^), the polarizability term (see 
Equation |6]l should include a quadrupole term in addition 



o 
o 




1.0 lO'' 
Energy (eV/u) 

Fig. 2. Total (thick solid curve) and final-state resolved cross sec- 
tions for CT recombination of Kr^^ . Energies are given in eV/u, 
where u is the atomic mass unit. 



to the dipole polariz ability if X is not an S ground term 
(lGentrv&Gieselll977h : 



ad 



ad 

2/;4 



_Q_ 

2R^' 



(11) 



To determine AH(R) for each element, static average elec- 
tric dipole polarizabilities were taken from iLid T (l200l 
and quadrupole mat rix elements from Tables 2 and 3 of 
iFroese Fisched (Il977l) . 

The transit ion prob ability for the reaction is p - e " with 
(lDemkovlll96 4: Swartz 1994t) 



1 + exp 



2n^AU(Rx) 
hav 



-1 



(12) 



where h is Planck's constant and a is the polarizability of the 
neutral species X. Defining 



^ _ TT^AUjR,) 
hav 

the cross section at energy E is 



o-{E) = kR: 



4-e- 



(1 -t-e-*-')^ 



dx. 



(13) 



(14) 



We compute two-channel cross sections, which can be added 
to determine the total cross section. While this addition can lead 
to overestimated total cross sections (see Sect. 12.11 ). only Se and 
Br have more than one exoergic final channel, and their cross 
sections are typically small in the relevant temperature range. As 
for the MCLZ calculations, the cross sections are integrated over 
a Maxwellian velocity distribution to derive rate coefficients. 

3. Results 

In Table [T0 we provide information for the exit channels of our 
MCLZ and Demkov calculations, including energy defects AE, 



Table[T]is available in the online version of this article. 
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avoided crossing distances (or interaction distances, for singly- 
charged ions) R^, and the energy separation AU of the adia- 
batic curves at /?x- These parameters were used to compute the 
CT cross section for each final channel. In addition, we list the 
molecular symmetries common to the initial and final channels 
and whether the reaction is Type I or Type II. When a final state 
can connect to the initial state via more than one molecular sym- 
metry, the cross sections for each symmetry were added to de- 
termine the state-specific cross section. The state-resolved cross 
sections were added to determine the total CT cross section, and 
rate coefficients determined by integrating over a Maxwellian 
velocity distribution. 

The total CT recombination rate coefficients of each ion are 
given in Table |2] as functions of temperature, and are illustrated 
in Figure [3] We note that for ^ = 5, highly-excited Rydberg 
states can be favorable exit channels, but the reported level ener- 
gies in NIST for X"*^ ions are incomplete at those energies. Thus, 
CT into those highly-excited states could not be computed, and 
as a result the rate coefficients for q - 5 are lower limits. As seen 
in Fig. [3] the computed ^ = 5 rate coefficients are smaller than 
those for q = 4 with the exception of Ge^+. Since CT rate coeffi- 
cients generally increase with charge for ^ > 3, we recommend 
that q = 4 CT rate coefficients be used as rough estimates for Se, 
Br, Kr, Rb, and Xe g' > 5 ions. 

When the radial coupling CT rate coefficient from our 
MCLZ and Demkov calculations is sufficiently small, CT can 
proceed via other coupling mechanisms. Two such examples 
are radiative CT and spin-orbit (SO) coupling. For radiative 
CT, the canonical rate coefficient is taken to be 10 cm^^ s ' 
(lButIeretalJll977|)^ This value is uncertain by at least an or - 
der of magnitude feutler et al.ll977HStancil & Zvgelmanll996h . 
but given the small rate coefficient, CT is unUkely to be an im- 
portant recombination mechanism compared to radiative or di- 
electronic recombination when radiative CT dominates radial 
coupling. The rate coefficients listed in the tables are the max- 
imum of the radially-coupled and radiative CT rate coefficients 
(thus, all values of 1 .OOE- 14 indicate that radiative CT dominates 
radially-coupled CT). For some exit channels, radiative CT dom- 
inates radial coupling at all temperatures from 10^-10^ K. States 
with such small rate coefficients tend to be associated with small 
or large avoided crossing distances (7?x < 5, > 15-20 qq). 
For some singly- and doubly-charged ions, the only exothermic 
exit channel(s) has a very small rate coefficient dominated by 
radiative CT. 

Spin-orbit coupling is not as well understood as the radial or 
radiative mechanisms. Rate coefficients for SO-coupled CT re- 
actions can range from 10~'^-10~^^ cm^ s~' (Butler & D algarnd 
119791 11980a; Pra dhan & Dalgarnolll994l: [Stancilll2001l) . though 
detailed calculations are very few and to our knowledge exist 
only for singly-charged ions. Because of this, we consider SO 
coupling only for CT recombination of singly-charged ions and 
CT ionization of neutral species. For small radial coupling CT 
rate coefficients (< lO^'^cm'' s"'), we adopt 10"'^ cm-' s"', 
which is the largest compu ted rate coefficient for SO coupling 
(iPradhan & Dalgarnoll 19941) . as an upper Umit to the rate coef- 
ficient. This estimate provides a much smaller, but more robust 
upper limit than estimates based on the Langevin model (Stancil 
120011) . 

However, SO coupling is not relevant for all of the singly- 
ionized species we consider, since it only connects initial and 
final states of different molecular symmetries for which radial 
and rotational coupling are forbidden. For CT recombination of 
Rb^ (and the reverse reaction), the entrance and exit channels 
correspond to molecular states, and SO coupling therefore 
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Fig. 3. Total rate coefficients for CT recombination of X^ (dash- 
dot-dot-dot-dash curves), X^^ (dash-dot-dash curves), X^^ (dot- 
ted curves), X'*^ (solid curves), and X^^ (dashed curves), where 
X is the element indicated in each panel. The ^ = 5 rate coef- 
ficients are lower limits due to incomplete experimental energy 
level measurements. Note that the Rb^ rate coefficient is below 
10 ''^ cm^ s ' over the plotted temperature range. 



plays no role. In the case of Kr^ and Xe^, for which there is 
only one exothermic exit channel, the diabatic potentials for the 
triplet exit channels do not cross the potential curve of the sin- 
glet entrance channel, and hence SO coupling is unlikely to be 
important for these reactions. We therefore adopt the SO cou- 
pling upper limit to the rate coeflScients only for singly-ionized 
Ge, Se, and Br 

CT ionization is exothermic for the neutrals of the elements 
we consider, with the exception of Kr. We provide rate coeffi- 
cients for these systems in Table [3] and depict them in Figure |4] 
As for CT recombination, we adopt the largest of the rate coeffi- 
cients associated with the radial, radiative, and (when important) 
SO mechanisms. 

The CT recombination rate coefficients for singly-ionized 
species were computed from the CT ionization rate coefficients 
via detailed balance: 



Q'rec ~ ^'m 



/70,rec / A£, 

1 exp 

PO.ion \ 



(15) 



where aiecjon are the rate coefficients for CT recombination and 
ionization, respectively, po.i-ecion are the approach probabilities 
corresponding to the molecular symmetries involved in the re- 
combination and ionization reactions, A/ioo is the asymptotic 
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Table 3. Total charge transfer ionization rate coefficients (cm^ s ') for the reactions X° + — > + H + 6E. 













7-(K) 












l.OE+02 


2.0E+02 


4.0E+02 


6.0E+02 


8.0E+02 


1.0E-H03 


2.0E-I-03 


4.0E-I-03 


6.0E-I-03 


8.0E-I-03 


Projectile 


l.OE+04 


2.0E+04 


4.0E+04 


6.0E+04 


8.0E+04 


1.0E-H05 


2.0E+05 


4.0E-H05 


6.0E-h05 


8.0E-h05 


Ge" 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 




l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


1.18E-12 


2.32E-12 


5.83E-12 


Se" 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


l.OlE-12 


1.02E-12 




1.04E-12 


1.86E-12 


1.24E-11 


4.24E-11 


9.51E-11 


1.71E-10 


8.41E-10 


3.11E-09 


6.00E-09 


9.19E-09 


Br" 


l.OlE-12 


l.OlE-12 


l.OlE-12 


1.02E-12 


1.04E-12 


1.09E-12 


2.37E-12 


1.48E-11 


4.42E-11 


9.07E-11 




1.52E-10 


6.28E-10 


1.98E-09 


3.51E-09 


5.05E-09 


6.58E-09 


1.36E-08 


2.50E-08 


3.45E-08 


4.29E-08 


Ki« 


3.62E-45 


1.80E-31 


1.42E-23 


1.62E-20 


8.27E-19 


1.09E-17 


4.93E-15 


3.27E-13 


2.14E-12 


6.70E-12 




1.48E-11 


1.17E-10 


5.70E-10 


1.22E-09 


1.96E-09 


2.75E-09 


6.81E-09 


1.41E-08 


2.04E-08 


2.59E-08 


Rb° 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 




l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OlE-14 


1.12E-14 


Xe" 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


l.OOE-14 


1.04E-14 


1.48E-14 


3.34E-14 




8.32E-14 


1.60E-12 


1.91E-11 


6.5 IE- 11 


1.42E-10 


2.49E-10 


1.13E-09 


3.78E-09 


6.84E-09 


9.99E-09 



Notes. For Ge, Se, and Br, the spin-orbit coupling rate coefficient found for Cl^, 10"'- cm' s"' dPradhan & Dalgamolll994t) . is adopted as an upper 
limit when the rate from radial coupling falls below that value. SO coupling is not relevan t for CT rea ctions involving neutral Kr, Rb, and Xe, and 
hence we set the rate coefficient to the canonical radiative CT rate (10"'"* cm' s~': . Butler et al.lll977t) when the radial coupling rate coefficient is 
smaller than that value. Note that radiative CT does not occur for the endoergic reaction Kr" -I- — » Kr^ + H + 6E. 



in Table 4 and for CT ionization in Table 5 (both available at the 
CDS). 

The derived abundances of Se and Kr have explicitly been 
shown to be sensiti ve to uncertainties in CT rate coefficients 
(ISterling et alj|200"7l) . and other trans-iron elements are expected 
to be similarly sensitive. Unfortunately, it is not possible to as- 
certain robust uncertainties for CT rate coefficients computed 
with the LZ or Demko v approximations wi t hout ac companying 
QMOCC calculations. iButler & Dal^arnd (Il980bl) found that 
their LZ rate coefficients agreed with quantal calculations to 
within a factor of three for ions of light elements with large CT 
rate coefficients, but that the disagreement could exceed an order 
of magnitude in the case of systems with small rate coefficients. 
No such comparison is possible for the heavy ions we have in- 
vestigated, since QMOCC computations are not available for any 
of these systems. Our forthcoming study of the sensitivity of n- 
capture abundance determinations to atomic data uncertainties, 
using Monte Carlo simulations with photoionization modeling 
codes, will help to identify the ions that most urgently require 
QMOCC determinations of CT rate coefficients. 

4. Summary 

We have presented CT recombination rate coefficients for the 
first five ionization stages of the «-capture elements Ge, Se, 
Br, Kr, Rb, and Xe reacting with neutral hydrogen. A multi- 
channel Landau-Zener approach was employed for multiply- 
charged ions, while the Demkov approximation was used for 
singly-charged species. We note that the computed rate co- 
efficients for q = 5 are lower limits, due to the incom- 
pleteness of energy listings in the NIST database. Because of 
that, we recommend that the four-times ionized recombination 
rate coefficients be used for <7 > 5, with the exception of 
Ge^^. For reactions between the neutral species of these ele- 
ments and H^, we computed CT ionization rate coefficients. 
Total rate coefficients are given over the range of tempera- 
tures 10^-10^ K in Tables 12] and [3] while final-state resolved 
rate coefficients are presented in Tables 4 and 5 (available at 
the CDS). In addition, all rate coefficients can be obtained 
at http://www.pa.msu.edu/astro/atomicdata/charge_transfer/ In 
tandem with our concurrent investigations into the photoion- 
ization and radiative and dielectronic recombination properties 
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Fig. 4. Total rate coefficients for CT ionization of the indicated 
neutral species. 



energy defect, and the Boltzmann constant. We note that 
Equation [15] was applied only to state-resolved rate coefficients. 
In the case of Kr, for which CT recombination of singly-charged 
ions is exothermic, the CT ionization rate coefficients were com- 
puted in an analogous manner. For neutral Kr, the energy de- 
fect is very small (corresponding to 4700 K), and hence the 
CT ionization reaction can operate efficiently in photoionized 
nebulae. CT recombination of singly-charged ions is less im- 
portant at photoionized temperatures for the other elements (see 
Fig. [3]), as the lowest energy defect is that of Xe (corresponding 
to 17000 K). 

State-resolved rate coefficients are important for understand- 
ing the emission line spectra of photoionized nebulae, since 
CT populates e xcited states fro m which emission lines can be 
measured (e.g., iKingdon & Ferla nd 1996). We therefore pro- 
vide final-state resolved rate coefficients for CT recombination 
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of w-c apture elements (e.g.. lSterling & Witthoeftir201 lUSterlingl 
12011"), these CT rate coefficient determinations will enable the 
abundances of «-capture elements to be reliably determined in 
ionized astrophysical nebulae. Such studies bear significant im- 
plications for n-capture nucleosynthesis, the structure and inter- 
nal mixing of evolved stars, and the chemical evolution of trans- 
iron elements in the Universe. 
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Table 2. Total charge transfer recombination rate coefficients (cm^^ s"') for the 
reactions X''^ + H — » X*''"'*+ + + SE. Each listed rate coefficient is the largest 
of the radially-coupled CT rate coefficient from our MCLZ and De mkov calcu- 
lations and the canon ical radiative CT rate coefficient 10"''' cm^ s"' dButler et alj 
Il977l : lstancil & Zyge lman 1996). Note that radiative CT does not occur for en- 
doergic reactions. 
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l.OOE-14 


l.OOE- 14 


l.OOE- 14 


l.OOE-14 


2.00E- 


■14 


5.90E- 


-14 


1.37E- 


■13 




2.68E-13 


2.17E-12 


1.54E-11 


4.39E-11 


8.78E-11 


1.46E-10 


6.08E-10 


2.03E- 


■09 


3.74E- 


-09 


5.57E- 


■09 



7 
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Table 1. Details of MCLZ and Demkov calculations, including exit channels, molecular symmetries, energy defects A£, avoided crossing 
interaction distances R^, and the energy separations of the adiabatic potentials AU(Rx)- The molecular symmetries of each entrance channel 
given in parentheses after the reaction. 



Exit 


AE 




AU{R,) 


Molecular 




Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries'' 


Type 


Ge"(3p) + H+('S) 


^Ge^ 


- + H(-S) + AE 


{'I.-, 3n) 




4j./ 4p zp" 


0.2041 


7.258 


0.2089 


3s- 


I 


4s2 4p 2po 


0.2041 


7.406 


0.1943 


3n 






('S) + H(^S)^Ge+ +H+('S) + A£ (-1;+) 




4s^ 4p -P" 


0.0805 


12.48 


5.922E-03 




I 


Ge3+(-S) + H(-S) - 


->Ge^+ 


+ H+('S) + A£ 


('s+,3i;+) 




4s4p3p° 


0.4633 


4.648 


0.2070 


3l+ 


I 


4s 4p 'P" 


0.3333 


6.240 


7.592E-02 




I 


Ge-'+ 


('S) + H(-S)^Ge 


3++H+('S) + A£ (-Z+) 




4p^P° 


0.8013 


4.138 


0.2732 




I 




0.3111 


9.768 


5.463E-03 


2S+ 


I 


5s 2S 


0.2724 


11.11 


1.847E-03 


2S+ 


I 


Ge=+(^D) + H(-S) Ge'*+ 


+ H+('S) + A£' 


'n 'A 3z+ 


3n, 3 A) 


3d''4d3s 


0.7481 


5.723 


0.1071 


3i;+ 


I 


3d''4d 'P 


0.7219 


5.901 


9.530E-02 


'n 


I 


3d''4d3G 


0.7202 


5.913 


9.454E-02 


3i;+,3n,3^ 


I 


Sd' 4d 3p 


0.7161 


5.942 


9.274E-02 


3n 


I 


3d' 4d 3d 


0.7095 


5.990 


8.983E-02 


3i;+, 3n, 3a 


I 


3d'4d 'F 


0.6993 


6.065 


8.544E-02 


'n.'A 


I 


3d'4d 'D 


0.6905 


6.132 


8.168E-02 


'n, 'A 


I 


3d'' 4d 'F 


0.6903 


6.134 


8.157E-02 


3n, 3a 


I 


3d'' 5s 3d 


0.5888 


7.062 


4.274E-02 


3z+, 3n, 3a 


I 


3d''4d 'S 


0.5763 


7.200 


3.870E-02 


'1+ 


I 


3d'' 5s 3p 


0.5677 


7.299 


3.603E-02 


3n 


I 


3d'' 5p 3p'' 


0.4370 


9.313 


7.827E-03 


3z+,3n 


I 


3d'' 5p 3F'' 


0.4313 


9.430 


7.139E-03 


3i;+, 3n, 3a 


I 


3d' 5p 'D" 


0.4277 


9.506 


6.724E-03 


'n.'A 


I 


3d' 5p 3d'' 


0.4149 


9.786 


5.385E-03 


3n, 3a 


I 


3d' 5p 'P" 


0.4136 


9.815 


5.262E-03 


'n 


I 


3d'5p 'F" 


0.4107 


9.882 


4.989E-03 


'n, 'A 


I 


Exit 


AE 


Rx 


Af/(Rx) 


Molecular 




Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries'' 


Type 


Se"(3p) + H+('S) 


^Se+ 


+ H(^S) + AE 


(3^-, 3n) 




4s^ 4p3 "S" 


0.1413 


7.006 


0.1304 


3S- 


I 


4s2 4p3 ^D" 
4s2 4p3 ^D" 
4^,2 4p3 2po 


0.0796 


7.860 


7.958E-02 


3n 


I 


0.0796 


8.020 


7.240E-02 


3s- 


I 


0.0338 


9.251 


3.396E-02 


3n 


I 


Se^+(3p) + H(-S) ^ Se+ + H+('S) + AE (^S", ^n, "S", *' 


n) 


4j,2 4pJ 2po 


0.1714 


6.040 


8.688E-02 




I 


Sgj+(2pO) ^ 


H(-S) Se-^* + H+('S) + AE 'n, 3^+, 


3n) 


4s^4p^ 'S 


0.5034 


4.256 


0.2568 


's+ 


I 


4s4p3 3d° 


0.2037 


9.921 


4.588E-05 


3n 


II 


Se-'+ 


('S) + H(^S)-^Se 


'+ + H+('S) + A£ 




4s 4p'^ -^D 


0.6023 


5.328 


8.177E-03 


^s+ 


II 


4s4p2 2s 


0.4917 


6.374 


2.706E-03 


2S+ 


II 


4s2 4d ^D 


0.3793 


8.088 


2.010E-02 


2S+ 


I 


4s2 5s 2 S 


0.3620 


8.452 


1.525E-02 


2S+ 


I 


4s2 5p ^P" 


0.2095 


14.38 


1.178E-04 


2S+ 


I 


Se'+(-S) + H(-S) - 


Se''+ 


+ H+('S) + A£ 


('S+,3S+) 




4s4d 'D 


1.0394 


4.342 


0.2453 


is+ 


I 


4s4d3D 


0.8359 


5.209 


0.1483 


3S+ 


I 


4s4d3D 


0.7012 


6.051 


8.625E-02 


3S+ 


I 


4s4f 'F" 


0.3639 


11.11 


1.854E-03 


's+ 


I 


Exit 


AE 


Rx 


A(7(Rx) 


Molecular 




Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries'' 


Type 


Br"(^P") + H+('S)T±Br 


^ + H(^S) + AE 


(^s+, ^n) 




4sMp'' 3p 


0.0589 


7.241 


6.159E-02 


% 


I 


4s2 4p'' 'D 
4s2 4p'' 'D 


0.0105 


9.612 


1.191E-02 




I 


0.0105 


10.05 


8.681E-03 


2S+ 


I 


Br^+C*S'') + H(^S) 


-*Br+ 


+ H+('S) + AE 


(3S-, ^S") 




4s^ 4p'' 3p 


0.2947 


3.576 


0.3579 


3s- 


I 


Br3+(3p) + H(^S) ^ Br^+ + H+(' S) + A£ (^S", ^n, '^'S,-, 




4s Ap" ^D 


0.4068 


5.219 


9.140E-03 




II 


4s2 4p2 5s "P 


0.1462 


13.73 


2.046E-04 


*s-, "n 


I 


4s2 4p2 5s 2p 


0.1346 


14.91 


7.473E-05 


2s-, 2n 


I 
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Table 1. Continued. 







Br'+ + H+('S) + A£ 


fly 

\ z. 


+ iTT 3y+ 
, il, z. , 






U.J /Zo 


5.570 


6.369E-03 




11 




Ac- A n3 3 DO 

4s 4p r 


0.5091 


6.176 


3.352E-03 




3v+ 3rr 




4s4p3 ^D'^ 


U. jUUz 


6.277 


3.006h-03 




11 




4S 4p 4a r 


n /IT 1 T 


7.336 


3.5t)7h-02 




3v+ 3tt 
h , 11 




Af A ri3 1 T)(3 

4S 4p r 


U.4i J / 


/.4Zo 


Q /IT 1 c r\A 
o.4ziii-U4 




1 V+ 1 TT 
2i , 11 








8.074 


2.031E-02 




3^ 




4c2 4n4d 3po 


0.3468 


8.803 


1.165E-02 










0.3448 


8.852 


1.121E-02 




ly+ Itt 
Zj , 11 






0.3301 


9.227 


8.373E-03 




ly+ Itt 
Zj , 11 




4s2 4p 4d ' D" 

/lt.2 /In 3po 
H-s Hp J& 1 


0.3245 


9.380 


7.426E-03 




'n 




0.3180 


9.564 


6.422E-03 








4s2 4p5s ipo 


0.3035 


10.04 


4.526E-03 




'n 




4s2 4p5p 


0.1619 


18.57 


2.983E-06 










('S) + H(^S) ^ Br''+ + H+('S) 4 


- A£ 






4s 4p'^ 


1.1321 


4.040 


2.853E-02 






jj — 


4s2 4d 


0.8360 


5.209 


0.1483 








4s2 5s2s 


0.7214 


5.904 


9.511E-02 








Exit 


A£ 


Rx 


AJ/(Rx) 


Molecular 




Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries'' 


Type 


Kr+C^P") + HC^S) Kr + H+C'S) + A£ 




'n, ^Y.^, ^' 


n) 


4s^ 4p*' ' S 


0.0147 


8.394 


1.595E-02 




Li 




Kr^+C""?) + H(^S) -> 


Kr+ +H+('S) + A£ 




Ztt 4v- 
, 11, L , 




4s2 4p5 ip> 


0.3873 


2.677 


0.4928 




11 






+ H+('S) + ^E 




— 


48^?"^? 


0.8485 


2.577 


0.5047 








4s2 4p3 4d5D" 


0.2269 


8.940 


1.047E-02 








4s2 4p3 5s5S° 


0.1942 


10.39 


3.313E-03 








4s2 4p3 4d ^D" 


0.1816 


11.10 


1.870E-03 








4s2 4p3 5s3s° 


0.1675 


12.01 


8.767E-04 








4s2 4p3 4d3G" 


0.1264 


15.86 


3.248E-05 




35;- 




4s2 4p3 5s3D'' 


0.1099 


18.23 


4.025E-06 




3y- 
z 




Kr^+C^P) + H(-S) -» 


Kr^+ +H+('S) + AE 




- irr 4v- 
, tt, z , 


it) 


4s 4p'' "P 


0.8787 


3.834 


3.429E-02 




L , 11 




4s2 4p2 4d -P„ 


0.6847 


4.775 


0.1924 




Z. , 11 




4s2 4p2 4d''F 


0.6265 


5.155 


0.1533 




4^— 4tt 
\ , 11 




4s2 4p2 4d''D 


0.6030 


5.328 


0.1378 




4tt 
11 




4s2 4p2 4d 


0.5937 


5.400 


0.1317 




2v- 2rT 
L , 11 




4s 4?-* 2p 


0.5532 


5.744 


5.308E-03 




2v- 2rT 
L , 11 




4s2 4p2 4d'*P 


0.5079 


6.193 


7.838E-02 




4v— 4rT 
L , 11 




4s2 4p2 5s "P 
4s2 4p2 4d2D„ 


0.4901 


6.394 


6.834E-02 




4^— 4rT 

2; , 11 




0.4719 


6.615 


5.864E-02 




2rr 
11 




4s2 4p2 5s 2p 


0.4691 


6.650 


5.723E-02 




2y- 2tt 

Z/ , 11 




4s2 4p2 4d ^Dfc 


0.4386 


7.069 


4.253E-02 




2n 




4s2 4p2 4d ^Pi 


0.4257 


7.266 


3.690E-02 








4s2 4p2 4d ^4 


0.4240 


7.293 


3.618E-02 








4s2 4p2 5s 


0.4124 


7.482 


3.152E-02 




2n 




4s2 4p2 4d -Dr 


0.3710 


8.259 


1.766E-02 




2n 




4s2 4p2 5p-S'' 


0.3595 


8.508 


1.461E-02 








4s2 4p2 5p ''D" 


0.3289 


9.260 


8.159E-03 








4s2 4p2 5p ■♦P" 


0.3203 


9.498 


6.766E-03 




'^n 




4s2 4p2 5p ■♦S" 


0.3058 


9.931 


4.797E-03 








4s2 4p2 5p ^D" 


0.3036 


10.00 


4.540E-03 








4s2 4p2 5p 2p° 


0.2884 


10.51 


3.011E-03 




2n 




4s2 4p2 5p ^F" 


0.2543 


11.88 


9.761E-04 




2n 




4s2 4p2 5p 2d° 


0.2512 


12.03 


8.658E-04 




21;-, 2n 




4s2 4p2 5p 2p° 


0.2235 


13.49 


2.522E-04 




2n 




Kr^+C^P") + 


H("S) 


Kr''+ +H+('S) + A£ 


('X 


+, 'n,^'i+. 


it) 


4s^4p4d 'P" 


0.9083 


4.841 


0.1851 








4s2 4p4d3D" 
4s2 4p4d 'F° 


0.8781 


4.998 


0.1687 




3n 




0.8110 


5.344 


0.1364 




'1;+, 'n 




4s2 4p4d 'P° 


0.7946 


5.437 


0.1287 




'1;+, 'n 




4s2 4p5s3p'' 


0.7633 


5.619 


0.1146 








4s2 4p5s 'P° 


0.7342 


5.816 


0.1008 




'1;+, 'n 




4s2 4p5p 'P 


0.5858 


7.094 


4.177E-02 




'n 




4s2 4p5p 


0.5761 


7.202 


3.865E-02 








4s2 4p5p 3p 


0.5554 


7.446 


3.237E-02 




3n 




4s2 4p5p3S 


0.5395 


7.646 


2.794E-02 








4s2 4p5p 'D 


0.5264 


7.822 


2.452E-02 




'1;+, 'n 




4s2 4p5p 'S 


0.4759 


8.592 


1.370E-02 








Exit 


A£ 




Af/(Rx) 


Molecular 





N 
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Table 1. Continued. 



Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries" 


TvDe 


Rb"(^S) + H+('S)T±Rb^ 


- + H(^S) + 


AE 


^ ) 




4s'^ 4p'' ' S 


0.3462 


10.46 


0.3332 




2y + 


J 


Rb^+(-P'') + H(^S)- 


>Rb+ +H+('S) + A£ 


( ^ 


f iTT + 




4s^ 4p*' ' S 


0.5031 


2.175 


0.5374 




lv+ 
2- 




Rb^+('P) + H(-S) ^ 


Rb-+ +H+('S) + A£ 


(5; 




it) 


4s^ Ap" 4d ''D 


0.2275 


8.917 


1.066E-02 




4n 
11 




4s2 4p* 4d ''F 


0.1798 


11.20 


1.710E-03 






J 


4s2 4p'* 5s 2d 


0.1629 


12.34 


6.637E-04 








4s2 4p'* 4d 2p 


0.1593 


12.62 


5.269E-04 








4s2 4p'' 4d 


0.1559 


12.89 


4.192E-04 








4s2 4p'' 4d 


0.1463 


13.73 


2.062E-04 








4s2 4p'' 4d 


0.1428 


14.06 


1.552E-04 








4s2 4p'* 5s 2p 


0.1261 


15.90 


3.142E-05 








4s2 4p'* 4d 


0.1219 


16.45 


1.950E-05 




2n 





Rb^'+C^'S") + H(-S) 


Rb3+ 


+ H+('S) + A£ 




) 


4s^4pJ 4d^D" 


0.6392 


5.043 


0.1642 








4s2 4p3 4d 3d;; 


0.5865 


5.458 


0.1270 








4s2 4p3 4d 3po 
4s2 4p3 4d3G° 


0.5553 


5.724 


0.1070 








0.5237 


6.027 


8.764E-02 








4s2 4p3 5s^S° 


0.4923 


6.364 


6.976E-02 




5s- 




4s2 4p3 5s3s° 
4s^ 4p3 4d ^D" 


0.4619 


6.744 


5.357E-02 




3s- 




0.4225 


7.316 


3.558E-02 




3s- 




4s2 4p3 5s3D° 


0.3955 


7.778 


2.534E-02 




3S- 




4s2 4p3 4d3s'' 


0.3935 


7.815 


2.465E-02 




3S- 




4s2 4p3 4d 3d° 
4s2 4p3 5p 5p 


0.3931 


7.822 


2.452E-02 




3S- 




0.3123 


9.732 


5.622E-03 




5S~ 




4s2 4p3 5p 3p„ 


0.2901 


10.45 


3.161E-03 




3S~ 




4s2 4p3 5p 
4s2 4p3 5p 3p,, 


0.2189 


13.77 


1.988E-04 




3S~ 




0.1942 


15.50 


4.465E-05 




3S- 




Rb'+(3p) + H(-S) ^ Rb''+ + H+('S) + AE 




'n, *s- 






1.3832 


3.388 


4.999E-02 








4s 4p'' 2d 


1.2509 


3.706 


3.833E-02 




2n 




4s 4p'' 2p 


1.1465 


3.996 


2.969E-02 




2s-, 2n 




4s2 4p2 5s "P 


0.7897 


5.465 


0.1264 




*s-, 




4s2 4p2 5s 2p 


0.7817 


5.513 


0.1226 




2s-, 2n 




4s2 4p2 5s ^D 


0.7110 


5.979 


9.049E-02 




2n 




Exit 


AE 


Rx 


AJ/(Rx) 


Molecular 




Channels" 


(a.u.) 


(a.u.) 


(a.u.) 


Symmetries 


" Type 


Xe"('S) + H+('S) ^ Xe+ + H(^S} + A£ 


('S+) 




5s2 5p^ 2po 


0.0380 


7.933 


3.515E-02 




is^ 




Xe2+(3p) + H(-S) ^ Xe+ + H+('S) + AE 




.'n.^s- 




5s2 5p^ 2p, 


0.2551 


4.078 


0.2817 




2s-, 2n 




Xe^+C'S") + H(-S) 


Xe2+ 


+ H+('S) + A£ 


(3s-, 'S- 


) 


5s2 5p3 5d ^D" 


0.1307 


15.35 


5.097E-05 




is- 




Xe*+(3p) + H(-S) -» Xe3+ + H+('S) + AE 




,2n, *s- 





5s5p^^P 


0.5315 


5.947 


4.283E-03 


^s-, ''n 11 


5s 5p'' 2d 


0.4386 


7.068 


1.259E-03 


2n 11 


5s2 5p2 5d 2p„ 


0.3920 


7.843 


2.414E-02 


2s-, 2n 1 


5s2 5p2 5d "F 


0.3608 


8.479 


1.494E-02 


^s-, "n 1 


5s2 5p2 5d 2f„ 


0.3494 


8.741 


1.222E-02 


2s-, 2n 1 


5s2 5p2 5d "tD 


0.3172 


9.587 


6.306E-03 


''n 1 


5s2 5p2 5d "^P 


0.2713 


11.15 


1.781E-03 


''s-,'*n 1 


5s2 5p2 5d 2g 


0.2647 


11.43 


1.425E-03 


2n 1 


5s2 5p2 6s "P 


0.2450 


12.32 


6.749E-04 


''s-,'*n 1 


5s2 5p2 5d 2d„ 


0.2437 


12.39 


6.397E-04 


2n 1 


5s2 5p2 6s 2p 


0.2237 


13.48 


2.548E-04 


2s-, 2n 1 


5s2 5p2 5d 2F4 


0.1969 


15.29 


5.361E-05 


2s-, 2n 1 


5s2 5p2 5d 2D6 


0.1935 


15.55 


4.253E-05 


2n 1 


5s2 5p2 5d 2p4 


0.1867 


16.12 


2.605E-05 


2s-, 2n 1 


5s2 5p2 4f''G'' 


0.1645 


18.27 


3.870E-06 


'^s-, ''n 1 


Xe^+(2p") H 


- H(2S) ^ 


Xe'*^ +H+('S) + A£ 


('s+, 'n, 3s+,3n) 


5s2 5p4f 3d 


0.9346 


4.745 


0.1958 


3s+, 3n 1 


5s2 5p4f 'D 


0.8902 


4.942 


0.1744 


is+, 'n 1 


5s2 5p6p3D 


0.7893 


5.468 


0.1262 


3s+, 3n 1 


5s2 5p6p3p 


0.7835 


5.502 


0.1235 


3n 1 


5s2 5p6p 'P 


0.7696 


5.586 


0.1170 


'n 1 


5s2 5p6p3s 


0.7487 


5.719 


0.1074 


3S+ 1 


5s2 5p6p 'D 


0.7361 


5.803 


0.1017 


is+, 'n 1 


5s2 5p6p 'S 


0.6947 


6.100 


8.346E-02 


lS+ 1 
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Table 1. Continued. 



" Exit channels refer to states in the resulting ion after CT with H or H^. Only information about electrons in the valence principal quantum 
number shell is given. Some terms are given a lettered sub-index (e.g., -D„) if multiple terms with the same configuration and orbital and spin 
angular momenta are exoergic exit channels. Note that exit channels with rate coefficients smaller than 10"'* cm^ s"' (the canonical radiative CT 
rate) in the temperature range 10^-10* K are not listed, unless no exoergic exit channels with larger rate coefficients exist. 
Only molecular symmetries that couple to the entrance channel are listed. 



